// Solve the Momentum equation

MRF.correctBoundaryVelocity(U);

fvVectorMatrix UEqn
(
    fvm::ddt(U) + fvm::div(phi, U)
  + MRF.DDt(U)  //mrf.ddt member:ddt,class mrf
  + turbulence->divDevReff(U)//turbulence->divDevReff(U)即扩散项，因为湍流模型，流变模型等代码实际上是对扩散项进行修正，因此在  OpenFOAM中扩散项的计算归于湍流模型的代码中
 ==
    fvOptions(U)
);

UEqn.relax();

fvOptions.constrain(UEqn);

if (piso.momentumPredictor())
{
    solve(UEqn == -fvc::grad(p));

    fvOptions.correct(U);
}
